Stochastic incoherence in the response of rebound bursters 
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At an optimal value of the noise intensity, the maximum variability in rebound burst durations 
is observed and referred to as a response stochastic incoherence. A general mechanism underlying 
this phenomenon is given, being different from those reported so far in excitable systems, ft is 
shown to be determined by (i) the monotonous reduction of the hysteresis responsible for bursting 
caused by noise and consequent transformation of responses from rebound bursts to single spikes, 
and (ii) a symmetry breaking in distributions of burst durations caused by existence of the minimum 
response length. The phenomenon is studied numerically in a Morris-Lecar model for neurons and 
its mechanism is explained with the use of canonical models describing hard excitation states. 
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PACS numbers: PACS: 05.45.-a, 05.10.-a 



INTRODUCTION 

o 

The presence of noise in nonlinear systems has been shown to induce non-trivial behavior in their dynamics 
Stochastic resonance Q occurs when the system is driven by both a periodic signal and noisy fluctuations. At 
an optimal value of the noise amplitude, the system exhibits the maximum correlation with the periodic signal. 
Oh. This phenomenon has been exhaustively studied theoretically and demonstrated experimentally in different kinds of 
nonlinear systems 0], e.g. in excitable j4| and bistable [B[ ones. It may occur, however, that even in the absence of 
periodic forcing the system reveals coherent oscillations at an optimum noise amplitude. This is the case of stochastic 
coherence (or coherence resonance) The reason for occurrence of stochastic coherence in excitable systems 

has been attributed to the existence of two different characteristic times |{|, activation and excursion times, which 
• <-h ' induce a coherence maximum in the inter spike intervals when they both reach a mutual minimum. Another proposed 
mechanism given in [Io| associates the appearance of coherence maximum with the noise induced bifurcation from the 
excitable to oscillatory regime. Stochastic coherence has been demonstrated theoretically in many systems including 
Q-ii excitable [l~D] and autonomous bursting systems flil ITii ] . It has also been demonstrated experimentally in various 
physical systems 

H ■ While stochastic resonance can appear in systems exhibiting various dynamical regimes, stochastic coherence can 
i> appear only in excitable systems, since it requires the existence of a refractory period. Excitability is a crucial feature 
; of neurons which enables an easy and efficient interaction between them [17] . In two-dimensional excitable systems, 
an above-threshold perturbation applied at time t = to, triggers a large excursion in phase space, until it finally 
returns back to the stable fixed point. During the spike duration the system is unsensitive to external perturbations. 
This temporal period marks a refractory period of the system. Higher-dimensional systems can give rise to a different 
type of excitability, namely, excitable or rebound bursting. The system being in the state of excitable bursting (at 
t > to), with the difference to standard excitable systems (e. g. FitzHugh-Nagumo model 

0). 

is susceptible to 

noisy fluctuations. Thus it has rather a pseudo-refractory period (pseudo, since it can be modified by the presence of 
noise). On the other side, in the absence of noisy fluctuations at t > to, the excitable burst has a well defined width 
determined by system parameters. 

In addition to the stochastic coherence, the phenomenon of stochastic incoherence (SI) has been reported, showing 
that in excitable systems exists a range of noise amplitudes at which the maximum inter-spike variability is observed. 
The phenomenon has been demonstrated in leaky integrate-and-fire model [l9| where the mechanism of occurrence 
has been attributed to the third time scale existing in excitable systems, the absolute refractory period. SI has been 
also shown to occur in FitzHugh-Nagumo model [20| where the underlying mechanism has been associated with the 
nonmonotonous relaxational behaviour of the system near the oscillatory regime. 

In this work a different type of stochastic incoherence is reported, namely, response stochastic incoherence (RSI). Its 
occurrence is found in a dynamical regime of excitable bursting. In the case of two dimensional excitable systems the 
spike duration (refractory period) is approximately fixed while the inter-spike interval is extremely sensitive to noise. 
In the case of bursting systems subject to noise, maximal variability is observed both in the inter-burst intervals and 
in burst durations. However, the two phenomena are observed in different range of noise amplitudes. In particular, 
the maximal variability in the inter-burst interval occurs for higher noise levels and refers to stochastic incoherence. 
In this study, the generation of bursts is triggered by a stimulus and not by the accompanying noise, whose only 
role is that of inducing variability in the duration of the elicited burst. It is shown, that the variability of excitable 



2 



burst durations Tg manifests a pronounced maximum at the optimal amplitude of fluctuations. The phenomenon is 
referred to as an incoherent in order to maintain a close relation with a well-known SI observed in excitable systems. 
The meaning of an incoherence in the case of RSI is associated with a maximum irregularity observed in Tg as the 
noise intensity is varied. It is demonstrated that the mechanism of the variability maximization in response duration 
is mainly determined by the sensitivity of the system hysteresis, which determines bursting regimes, to noise. During 
a decrease of the mean burst duration, or equivalently, a decrease of the hysteresis depth, the system transforms from 
hard excitation state (hysteretic) into a soft excitation state (non-hysteretic), leading to a single spike mean responses 
at larger noise amplitudes. The maximum variability occurs when the symmetry breaking in the distributions for 
burst durations takes place. This symmetry breaking is caused by the existence of the minimum possible response 
length (single spike) which becomes highly probable at larger noise amplitudes. The occurrence of the phenomenon is 
demonstrated in the Morris-Lecar model equations. Its underlying mechanism is explained on examples of canonical 
systems containing hysteresis. In particular, in models describing stabilized subcritical pitchfork bifurcation (or 
subcritical Hopf bifurcation), and in a model describing a damped pendulum. 



MODEL 

The general model equations are considered: 

{eix, e 2 y, e 3 z} = {f(x, y, z) + S{t), g(x, y), h(x, z)} (1) 

where x — y is the fast subsystem and z is the slow variable. The fast subsystem x — y gives autonomous spiking, 
meanwhile variable z rebounds spikes and enables bursting regimes, including excitable bursts. S(t) is an external 
perturbation acting as an external stimulus provoking the excitable response of a neuron and has the form of the Dirac 
S function. One chosen system parameter is then assumed to fluctuate following an independent white noise process 
with correlations (£(t)£(t')) = 2D6(t — t') and zero mean. The choice of the random process type is motivated 
by results reported in Ref. [21| which show that Gaussian noise has the same effect as channel or synaptic noise. 



Rebound bursting in Morris-Lecar system 

A special case of the general model in Eq. (JTJ is the modified Morris-Lecar model [22[ : 

etx = r(x) -2y(x + 0.7) -0.0Lz(x + 0.885) + S(t) 
e 2 y = C (- y + l/(l + e-<*-°-W 7 ) 

e 3 i = -z + 9.3(.t + 0.6) (2) 

where T(x) = -(1 W- ( * +0 - 01)/0 - 075) ) _1 (x-l)-0.065(x + 0.5), e 1 = 0.97, e 2 = l/cosh((x-0.1)/0.14) and e 3 = 30550. 
The Morris-Lecar model is a conductance-based model of voltage oscillations in the barnacle giant muscle fiber which 
is often used as a qualitatively accurate model of neuronal spiking. In order to obtain a rebound burst, the additional 
variable z is added to the subsystem x — y. An example of such a bursting response is shown in Fig. [T] (a). The 
system in Eq. ([5]) contains a dynamical hysteresis, i.e. a regime with transient bistability between a fixed point and 
repetitive spiking. The existence of such a dynamical hysteresis is responsible for bursting dynamics. The rest state 



disappears via fold bifurcation, and the periodic spiking disappears via saddle-homoclinic orbit bifurcation 23] . The 
latter bifurcation allows class I neural excitability, since it manifests in low frequencies appearing at the transition to 
a rest state onset. 

As already mentioned in the Introduction, the rebound bursting system has a rather pseudo-refractory period, since 
it can be influenced by the slight external fluctuations. On the other hand, the external stimulation has to be strong 
enough to elicit the rebound bursting response. Let consider the noiseless case when the system is stimulated with a 
rectangular single pulse of amplitude A and duration At. Numerical simulations show that the system responses have 
a form of a fix duration burst or single spike. The nonlinear responses characteristics to such external stimulation are 
shown in Fig. [1] (b). It can be seen that the response of the system is bistable (coexistence of single spike and fixed 
duration bursting responses) . A lack of variability in the burst durations and the two state responses of the system to 
stimuli with different parameters (A and At) suggest, that the system is rather unsensitive to fluctuations in external 
stimulation eliciting excitable burst. As a consequence, only fluctuations present in an upper spiking state can give 
rise to a variability in the response durations. In the presence of noiseless stimulus S, the system fires a burst of a 
fixed duration, with a well-defined passage from the lower to the upper state and vice versa, determined by internal 
system parameters. In Eq. @ the burst duration can be controlled by modulation of the parameter c. 
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RSI in rebound burster 



Two cases are considered: i) the parameter c is chosen to fluctuate c — > c + D£(t), modelling the internal noise 
or ii) the stimulus S is chosen to be accompanied by external fluctuations S — > S + D£(t). The rebound bursting 
can be measured in different manners. One way is to estimate the total time during which the repetitive spiking 
takes place. Another way is to count the number of spikes in a given burst, as usually done in experiments [24| . In 
the present work, the former way for measuring the burst durations Tg is used, namely, the estimation of temporal 
intervals during which the system fires spikes repetitively. In order to describe qualitatively the variability of the 
response durations, standard deviation is considered ot b = ((T B ) — {Tb) 2 )" 1 where () stands for mean. The statistical 
measurements are done on an ensemble of excited bursts in the presence of different realizations of noise. The noise 
intensities taken into consideration are small enough to keep FAR = (FAR- false alarm rate) at S — 0. This is done 
in order to study the specific response of the system to external stimulation. It is observed, that the mean burst 
duration decreases with increasing of noise intensity. In Fig. [2] (a) the effect of noise on the mean burst duration for 
different values of the parameter c is shown. This implies that in the presence of noise, the passage from the upper to 
the lower state becomes faster, and the transition time decreases as noise amplitude increases. The same behaviour is 
observed for both types of fluctuations present in the system. In view of the above results, it can be hypothesized that 
noise destroys the transiently existing bistability region in the bursting system, and transforms the hard excitation 
system into a standard excitable one-spike system or soft excitation system. This is clearly shown in Fig. [5] (b) where 
the distributions of the burst durations are plotted for different values of noise amplitude. From these distributions 
(Jt b is calculated and plotted in Fig. [3] (a) for various noise intensities D. The maximum of <7t b corresponds to a 
maximum variation in excitable burst durations. The choice of gt b to describe the variability of Tb instead of the 
usually used coefficient of variation CV = cf Tb / < T B > is caused by the fact that <T B > decreases with the increase 
of D. Thus, having a division by decreasing quantity, CV hides the useful information about Tb variability. In fact, 
it has been already pointed out in 25|, that CV sometimes fails as an indicator of coherence. 

It has been shown that the burst duration can be modified either by modulation of a control parameter c, either 
by varying the noise intensity. In fact, both, the effects of a control parameter c and noise amplitude D on the burst 
duration, are shown to be equivalent (see Fig. [3](b)). This implies that random fluctuations in the rebound burster 
can play a role of a control parameter. 



MECHANISM OF RSI 
Sensitivity of hard excitation states to noise 



In order to understand the effects of noise on the burst duration in the Morris-Lecar system let consider a canonical 
model of a system containing hysteresis which describes a stabilized subcritical pitchfork bifurcation [26| : 

x = rx + x 3 - x 5 (3) 

In the case, when the additional subsystem equation <j> = 2ir is added to Eq. ([3]) and x is considered as a radial 
variable, the system undergoes subcritical Hopf bifurcation [27j |. In the absence of the stabilizing term x 5 the system 
in Eq. (|3|) undergoes a supercritical pitchfork bifurcation when one unstable fixed point bifurcates to three fixed points, 
one stable and two unstable. The stabilizing term allows the re-occurrence of the two stable and one unstable fixed 
points. The common feature between the saddle- node (occurring in the Morris-Lecar model) and stabilized subcritical 
pitchfork bifurcations is the existence of low frequencies (in the case of saddle-node bifurcation) or low amplitudes (in 
the case of pitchfork bifurcation) appearing at the bifurcation onset. The two backward-bending branches of unstable 
fixed points bifurcate from the origin when r = tq = 0. Due to a stabilizing term x 5 , the unstable branches bend 
and become stable at r = r s , where r s < 0. The stable large-amplitude branches exist for all r > r s . In the range 
r s < r < 0, two different stable states coexist, the origin and the large-amplitude fixed points, which marks the 
hysteresis region. Inside this hysteretic region, the initial condition xo determines which fixed point is approached as 
t — y oo. In the white noise environment, the following stochastic differential equation is obtained for Eq. (|3|): 

dx t = (rx t + x 3 - x\ )dt + Dx t dW t (4) 

where the parameter r is assumed to fluctuate as r t — r + D^ t , where r is an average value, £t is a Gaussian white 
noise and D is the intensity of the fluctuations. Numerical simulations of Eq. (0| for different realization of noises 
reveal, that the same phenomenon as in the Morris-Lecar system occurs: the reduction of the mean hysteresis length 
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(corresponding to the burst duration in the Morris-Lecar system) with increasing noise intensity. The bifurcation 
diagrams for a system with and without noise are shown in Fig. [4] (a) and the distribution of the hysteresis lengths 
Ar for different noise amplitudes are plotted in Fig. [H (b). The mean hysteresis lengths (Ar) calculated from the 
distributions reveal the monotonous decrease until they reach (Ar) = 0. The mean values and its standard deviations 
a {Ar) arc shown in Fig. U (c). This implies the existence of a general mechanism of hysteresis reduction induced by 
noise. 

Another example of a system with hysteresis is the equation describing the dynamics of a pendulum with viscous 
damping (or a Josephson junction) [26j ] : 

4> = x 

x = fj, — sin (j) — ax (5) 

where a and [i are dimensionless control parameters. The parameter /i is assumed to fluctuate as [it = H + D£,t, where 
fi is an average value, £t is a Gaussian white noise and D is the intensity of the fluctuations. At small a and in the 
absence of noise (D = 0), as the parameter /i is increased the initially stable fixed point disappears in a saddle- node 
bifurcation at /i = 1 giving a limit cycle. If [i is brought back down, the limit cycle persists for \i < 1 and its frequency 
tends to zero as /i = \x c . This change in fi values as the bifurcation point is reached from different sides marks the 
hysteresis region. Now, when the internal parameter a is made fluctuating, the hysteretic region is seen to decrease 
with increasing noise amplitude (see Fig. 0] (d)). The same phenomenology of hysteresis reduction is observed also in 
this case (see Fig. H](e)). The averaged hysteresis depths (Afi) for different noise amplitudes are shown in Fig. H](f). 



Existence of a minimum mean response duration in the bursting system 

In the case of the Eqs. © and ([5]) the drift and diffusion processes continue until the averaged value (Ar) = and 
(A/i) = are reached, respectively. Then, when noise still increases, the drift disappears and only the amplitude of 
fluctuations around the states (Ar) = and (A/i) = increases (i.e. the diffusion increases with increasing D). In 
the case of a rebound burster, however, in which the hysteresis is a dynamic process allowing the transient spiking 
activity, the burst duration cannot be negative. Thus not only the drift stops but also the diffusion does. It occurs 
when the minimum possible burst duration T™ n , being a single-spike response like in a standard excitable system, 
is reached (see Fig. [5] (a)). The maximum variance is observed when the distribution is situated in the intermediate 
distance between two stable states: single spike and a maximum rebound burst defined at D = 0. Such a situation 
permits the system responses to visit all possible intermediate states giving the maximal variance. At larger noise 
amplitudes (but small enough to keep FAR — 0), the single spike responses become more probable giving rise to a 
Tb variability minimization. The value T B nm is a kind of a barrier which suppresses further increase of ot b ■ The 
variability increases with D until one tail of a distribution starts to disappear due to a collision with this barrier. In 
other words, the maximum variability of Tb is associated with the symmetry breaking in the form of distributions for 
Tg's consisting on the transitions from a normal to an exponential distribution form (see Fig. [2](b)). The symmetry 
of the probability distributions can be measured by the skewness coefficient 71, the standardized moment, defined as 
follows: 

71 = -3- (6) 

"Tb 

where ^ is the third moment about the mean (Tb) and <jt b is the standard deviation. The skewness of the normal 
distribution (or any perfectly symmetric distribution) is zero. The symmetry breaking in the distributions described 
by skewness coefficient 71 is shown in Fig. [5] (b). The correlation between minimum of 71 (tendency to a normal 
distribution form) and maximum of o~t b is clearly seen. 

In the case of Eq. [5] describing the dynamics of a pendulum with viscous damping, the distributions for the hysteresis 
depth tend to be symmetric as the noise amplitude increases. This is indicated by the skewness coefficient shown in 
Fig. |6l which at difference to the case of the bursting system, remains small and constant. 



DISCUSSION 



In this work the phenomenon of response stochastic incoherence in rebound bursters has been reported. The crucial 
ingredients for RSI to appear are the existence of hysteresis in the system and at least three dimensions, allowing a 
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fast-slow configuration and rebound bursting regime. The underlying mechanism has been provided, showing that the 
maximum variability in burst durations is caused by the sensitivity of the system hysteresis responsible for bursting 
dynamics to noise and by the symmetry breaking in distributions for the burst durations occurring as the barrier 
situated at the minimum value of the response duration (single spike response) is reached. The hysteresis vulnerability 
to noise has been demonstrated on examples of canonical systems describing hard excitation states. The implications 
of this finding is that the bursting responses to external stimulation are not necessarily well-determined and fixed, 
but depend strongly on the slight presence of random fluctuations. From the physical point of view the phenomenon 
of RSI here presented regards a different kind of stochastic incoherence, namely, the variance of the burst durations 
instead of interspike intervals as reported so far in the literature. 

The results on the hysteresis reduction could contribute to the understanding of some phenomena observed in 
models for bursting dynamics like the decreasing burst durations in the presence of noise [281 ] . On the other hand, 
RSI may be relevant in asynchronous firing patterns as already proposed for the case of SI in Ref. Q. 
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FIG. 1: (a) Excitable burst (upper time series) evoked by noiseless external stimulation S (lower time series plotted out of 
scale), (b) Response characteristics of the system to external stimulations of amplitude A and duration At. White color 
corresponds to a single spike response, and black color corresponds to fixed duration bursts. In (a) and (b) the Morris-Lecar 
system with c = 0.28 is considered. 
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FIG. 2: (a) Mean response durations for different values of the control parameter: c = 0.291 (solid line), c = 0.326 (dashed 
line) and c = 0.368 (dotted line), (b) Normalized distributions of the response durations in the case of c = 0.326 shown for 
increasing noise amplitudes (from right to left): C = {0,6 x 10~ 4 , 0.0014, 0.0022, 0.003, 0.0046, 0.0118} (the distributions for 
D = and D — 0.0118 are shown partially). 
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FIG. 3: (a) Standard deviation gt b versus noise amplitude D for different values of the control parameter: c = 0.291, c = 0.319, 
c = 0.347 and c = 0.368. (b) Dependence of (Is) on noise amplitude D and control parameter c. 
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FIG. 4: A canonical model for stabilized subcritical pitchfork bifurcation: (a) bifurcation diagram for a system without (dashed 
line) and with noise (solid line); the arrow marks the hysteresis reduction; (b) normalized distributions of the hysteresis lengths 
Ar for selected noise amplitudes (from right to left): D = {0.0474,0.1897,0.3795,0.5218,0.8064}; (c) mean hysteresis length 
(Ar) and its standard deviation cj / ai™) • Model describing damped pendulum: (d) response of the system to a variation of 
internal parameter \i (dashed line); upper trace corresponds to the system without noise and the lower to the system with 
noise; the arrow shows the hysteresis reduction; (e) normalized distributions of the hysteresis lengths A[i for selected noise 
amplitudes (from right to left): D = {0.0158,0.0632,0.1265,0.1897,0.3162}; (f) analogue to (c). 
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(a) 




FIG. 5: (a) Burst durations in the presence of noise: minimum possible T™ m (one spike) for D = 0.06, intermediate T B for 
D = 0.005 and in the absence of noise: T B determined by the value of parameter c, in this case c = 0.275. (b) Correlation 
between the standard deviation gt b (solid line, plotted out of scale) and skewness 71 (dashed line). Left panel corresponds to 
Eq. H with c = 0.284 and the right panel to c = 0.333. 




FIG. 6: Correlation between the standard deviation cta/j (solid line, plotted out of scale) and skewness 71 (dashed line) in the 
case of the equation describing the dynamics of a pendulum with viscous damping (Eq. [5J. 



